Self-Organization Induced Scale-Free Networks 
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What is the underlying mechanism leading to power-law degree distributions of many natural 
and artificial networks is still at issue. We consider that scale-free networks emerges from self- 
organizing process, and such a evolving model is introduced in this letter. At each time step, a new 
node is added to the network and connect to some existing nodes randomly, instead of "preferential 
attachment" introduced by Barabasi and Albert, and then the new node will connect with its 
neighbors' neighbors at a fixed probability, which is natural to collaboration networks and social 
networks of acquaintance or other relations between individuals. The simulation results show that 
those networks generated from our model are scale-free networks with satisfactorily large clustering 
coefficient. 
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The last few years have witnessed a tremendous activ- 
ity devoted to the characterization and understanding of 
complex networks 0, 0j j which arise in a vast number 
of natural and artificial systems, such as Internet t^l a l a , 
the World-Wide Web[|[|, social networks HHEmOI, 
airports network Il3lll4l. f ood webs [TH , biological in- 
teracting networks |l7lll8lll9l | and so forth. Particularly, 
much attention has been dedicated recently to the study 
of scale- free networks, namely, networks that display a 
power-law degree distributions, P(k) oc A; -7 , where k is 
the connectivity (degree) [2(j • And various proposals for 
dynamical evolution of scale-free networks has been intro- 
duced. Roughl y, t hese models can be classified into two 
main scenarios |2ll . One is under the mechanism "prefer- 
ential attachment" , which means new vertices are prefer- 
entially attached to existing vertices with large number 
of neighbors |U H3, 0, HU, and a related scenario is 
found in the protein duplication model j2|J. Another 
is a balance between a modelled tendency to form hubs 
against an entropy pressure towards a random networks 

However, what is the underlying mechanism leading to 
power-law degree distributions is still at issue. Differing 
from the two class of mechanisms mentioned above, we 
consider that self-organization maybe the fundamental 
mechanism which leads to power-law distribution of 
degree, and a model including such a mechanism is 
introduced in this letter. We think that the node newly 
added to the network will connect to some existing nodes 
randomly, not preferentially. For instance, considering 
the scientist collaboration networks, most of the times, a 
scholar comes into a new field not because he correlated 
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FIG. 1: A schematic representation of the evolving rules for 
the case m — 2, G — 2. At each time step, a new node is 
added to the network. For each new node, it will be active 
along G + 1 steps. In the first time step(i.e. the step when it 
was added to the network), it will randomly choose m nodes 
to be its neighbors. Then, in the following G steps, with 
probability p, this node will connect to the neighbors of the 
nodes which connected to it at the last step. To make it easier 
to understand, we draw the figure 1 as a sketch map for the 
case m = 2, G — 2. Assume node 7 connected to two nodes 
5 and (left) , then it will link to the neighbors of 5 and 6, 
i.e. 2, 3, 4i with probability p. Obviously, if some node is 
common neighbor of 5 and 6, the probability is 2p. At the 
same step a new node 8 is added to the network by linking 
to m~2 nodes randomly (omitted in Fig.l) (mid). Suppose 
node 7 connected to node 4 at last step, then it links to 
neighbors of 4, i- e - 3, 1, with the same probability p (right). 
The operation for each node lasts for G steps similarly. 



with some famous people in this field , but he read 
some experts' papers by chance and then connected 
with them or some of his partners were in this field 
already. Thus, a relatively new node having linked to 
some other nodes in the network will connect with its 
neighbors' neighbors, whose evidences can be found in 
collaboration and friendship networks. Therefore, we 
propose an evolving model following the rules: 
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FIG. 2: The statistical and evolutive characterization of de- 
gree, where m = 3,p = 0.06, G = 2, N = 2 x 10 4 : (a)degree 
distribution (b) normalized accumulative degree dis- 

tribution P(k) = kp(k)dk/k, (c)time evolution of max- 
imum degree fc max and (d)average degree fcavcragc- Obviously, 
p(k) oc k~"' , where 7 = 2.85 ± 0.07. Power fitness (a * t b ) of 
fcmax and fc aV erage are represented, with a = 5.06 ± 0.06, b = 
0.307 ± 0.005 and a = 5.42 ± 0.06, b = 0.018 ± 0.001 respec- 
tively, i.e. fc max oc TV 0,307 , fcavcragc oc jV - 018 . The number 
of edges of the networks E = JV * fc ^"°s° x ./V 1-018 , j n otner 
words, the networks are sparse. 
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1. Starting with a small number (too) of nodes which 
are global connected with each other. At each time step, 
a new node i is added to the networks, and connect to 
to nodes randomly. 

2. The node i will link with probability p to the 
neighbors of the nodes which connected to it at the last 
step. And at the next G — 1 steps (G is a fixed integer), 
let node i execute similar operation (as Fig.l shows). 

For simplification, the parameters to, p and G are con- 
stant. Obviously, rule*? is a self-organized process that 
impacts the topological structure of the networks. Fig. 2 
shows a typical scale-free network based on our model 
with N = 2 x 10 4 nodes. The exponent 7 is a little smaller 
than 3.0, which is close to empirical study. And time 
evolution of maximum degree /c max and average degree 
^average are exhibited in Fig. 2c and Fig. 2d, respectively. 
It can be found that fc max oc TV - 307 , fc avcragc oc N - 018 , 
thus the total number of edges in the network is pro- 
portionable to 7V 1+0 018 , in other words, the networks is 
sparse. Moreover, there is evidence to suggest that in 
real-world networks, e.g. World Wide Web, the average 
degree of nodes is increasing with time [3j. Our model 



FIG. 3: Degree distribution for different G and N: (a)G=l, 2, 
3, N = 10 4 ; (b) N = 10 4 ,2 x 10 4 ,3 x 10 4 , G=2, where m=2, 
p=0.06. It shows that the value of G and -/V impacts little the 
topology of network, i.e. different value of G and iV leads to 
almost the same power-law behavior of degree distribution. 



accords with that very well. 

In the following we will discuss how the parameters G 
and p impact the topology of networks. In fig. 3a one can 
see that for different G the model leads to almost the 
same degree distribution, i.e. the value of G impacts lit- 
tle the topology of networks. One can easily check that 
even for large G, which is close to real-world instances 
that people usually have wide communication with oth- 
ers, the networks based on our model are still of scale- 
free properties. This extensive result is very significant. 
However, when G is too large (e.g. G > 10) the com- 
puter is unable to work it out. Figure 3b exhibits that 
p{k) is independent of time and subsequently indepen- 
dent of network size N = niQ + t, which indicates that 
despite its continuous growth, the network self-organizes 
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FIG. 4: The degree distribution p(k) for networks with p=0.1, 
0.2, 0.3, where m=2, G=l, TV = 10 4 . It shows that for dif- 
ferent probability p, the degree distributions follows power- 
law behaviors p(k) oc k ' with almost the same exponent 
7 = 2.90 ± 0.09, and for larger p, the rang of power-law dis- 
tribution gets narrower. 
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FIG. 5: Clustering coefficient of networks with different size, 
TV, based on our model(triangle) and BA model(square). Ob- 
viously, the one of our model is much larger. 



to a scale-free stationary state. 

As illuminated above, for small values of p, e.g. 0.06, 
0.08 et al, the self-organization processes induce scale- 



free networks. However, when p is getting larger, that's 
to say, each node newly added to the network connects 
to its neighbors' neighbors with larger probability, the 
number of nodes with small degree will become fewer 
(as shown in Fig. 4). The results indicate that the rang 
of power-law degree distribution will get narrower as p 
increases. Let p=0, the network has exponential degree 
distribution. The analysis is simple. We can draw the 
master equation as follows: 

N{k,t+1) = 7V(fc,t)+ iV(fc ~ M) xm-^Mxm (1) 



t + mo 



t + mo 



where N(k,t) is the number of nodes with degree k at 
time step t. Since N(k, t) ~ p{k)t and t 3> mo for large 
t, seeking solution of this form we can get that p{k) oc 
e~m. How p impacts the degree distribution detailedly 
is significant and will be discussed in the future work. 

It is well known one of the shortcomings of BA's model 
is that the clustering coefficient is small and decreases 
with the increasing of network size, following approxi- 
mately a power law C ~ TV -0 75 lj. In Fig. 5, we report 
the clustering coefficient of the networks based on our 
model, which is much larger than that on BA's model. 
This is not difficult to explain. Suppose node hi and hi 
are neighbors of a common node /i3, then according to 
rule,2 mentioned above, node hi will link to node h2 with 
probability p if the former one is added to the network 
earlier than the later one. This induces that the edges 
among neighbors of node h3 will become more, thus the 
clustering coefficient will get larger. 

In summary, we propose a new model, based on "ran- 
domly attachment" and self-organization, which leads to 
a scale-free degree distribution of networks. The mecha- 
nism is different from preferential attachment mechanism 
where a scale- free distribution are generated during grad- 
ual growth of hubs. To investigate the "randomly attach- 
ment" mechanism is very important, since there are some 
networks that appear to have power-law degree distribu- 
tions, but for which preferential attachment is clearly not 
an appropriate modelflTl Il9l I30L l3l| . The simulation re- 
sults exhibit that, rather than BA networks, the networks 
generated from our model are of larger clustering, which 
is close to the real-life networks. Furthermore, we suggest 
that self-organizing process plays a major role in many 
real-life networks such as collaboration networks, social 
networks of acquaintance or other relations between indi- 
viduals and so on. We believe that many useful scale-free 
networks can be constructed using our approach. 



[1] R. Albert and A. -L. Barabasi, Rev. Mod. Phys. 74, 47 
(2002). 

[2] S. N. Dorogovtsev and J. F. F. Mendes, Adv. Phys. 51, 
1079 (2002). 

[3] M. E. J. Newman, SIAM Review 45, 167 (2003). 



[4] M. Faloutsos, P. Faloutsos and C. Faloutsos, Computer 

Communications Review 29, 251(1999). 
[5] R. Pastor-Satorras, A. Vazquez and A. Vespignani, Phys. 

Rev. Lett. 87, 258701(2001); R. Pastor-Satorras and A. 

Vespignani, Evolution and structure of the Internet: A 



4 



statistical physics approach (Cambridge University Press, [19] 
Cambridge, 2004). 

G. Caldarell, R. Marchetti and L. Pietronero, Europhys. [20] 
Lett. 52, 386(2000). 

R. Albert, H. Jeong and A. -L. Barabasi, Nature 401, 130 [21] 
(1999); A. -L. Barabasi, R. Albert and H. Jeong, Physica 
A281, 69(2000). [22] 
B. A. Huberman, The Laws of the Web (MIT Press, Cam- 
bridge, 2001); 

F. Liljeros, C. R. Edling, L. A. N. Amaral, H. E. Stanley 
and Y. Aberg, Nature 411, 907(2001). 

J. Scott, Social Network Analysis: A Handbook (Sage [23] 
Publications, London, 2000). 

S. Wasserman and K. Faust, Social Network Analysts 
(Cambridge University Press, Cambridge, 1994). [24] 
B. Hu, X. -Y. Jiang, J. -F. Ding Y. -B. Xie and B. -H. 
Wang, arXiv: |cond-m at/0408125 B. Hu, X. -Y. Jiang, 
B. -H. Wang, J. -F Ding, T. Zhou and Y. -B. Xie, arXiv: 
cond-mat/0408126 

A. Barrat, M. Barthlemy, R. Pastor-Satorras, and A. [25] 
Vespignani, Proc. Natl. Acad. Sci. USA 101, 3747(2004); 
A. Barrat, M. Barthlemy and A. Vespignani, Phys. Rev. [26] 
Lett., 92, 228701(2004).' 

R. Guimera, S. Mossa, A. Turtschi, and L. A. N. Amaral, [27] 
arXiv: cond- mat/0312535 

S. L. Pimm, Food Webs (University of Chicago [28] 
Press,Chicago,2002), 2nd ed. [29] 
A.E. Krause, K.A. Frank, D.M. Mason, R.E. Ulanowicz, 
and WW. Taylor, Nature 426, 282(2003). [30] 

H. Jeong, B. Tombor, R. Albert, Z. N. Oltvai and A. -L. 
Barabasi, Nature 407, 651(2000); H. Jeong, S. S. Mason, [31] 
A. -L. Barabasi, and Z. N. Oltvai, Nature411, 41(2001). 

J. Padani, Z. N. Oltvai, B. Tombor, A. -L. Barabasi and 
E. Szathmary, Nature Genetics 29, 54(2001). 



J. Stelling, S. Klamt, K. Bettenbrock, S. Schuster and E. 
D. Gilles, Nature 420, 190(2002). 

A. -L. Barabasi, Linked: The New Science of Networks, 
(Perseus Books Group, Cambridge MA, 2002). 

B. J. Kim, A. Trusina, P. Minnhagen, and K. Sneppen, 
arXiv: jcond-mat /0403006 

A. -L. Barabasi, and R. Albert, Science 286, 509 (1999); 
R. Albert and A. -L. Barabasi, Phys. Rev. Lett. 85, 
5234(2000); G. Bianconi and A. -L. Barabasi, Phys. Rev. 
Lett. 86, 5632(2001); G. Bianconi and A. -L. Barabasi, 
Europhys. Lett. 54, 436(2001). 

P. L. Krapivsky, S. Redner and F. Leyvraz, Phys. Rev. 
Lett. 85, 4629(2000); P. L. Krapivsky and S. Redner, 
Phys. Rev. E 63, 066123(2001). 

S. N. Dorogovtsev, J. F. F. Mendes and A. N. Samukhim, 
Phys. Rev. Lett. 85, 4633(2000); S. N. Dorogovtsev and 
J. F. F. Mendes, Europhys. Lett. 52, 33(2000); S. N. 
Dorogovtsev and J. F. F. Mendes, Phys. Rev. E 63, 
025101(2001). 

B. Tadic, Physica A293, 273(2001) 
A314, 278(2002). 

R. V. Sole, R. Pastor-Satorras, E. D. Smith, and T. Ke- 
pler, Advances in Complex Systems 5, 43 (2002). 
J. Berg and M. Lassig, Phys. Rev. Lett. 89 (22), 228701 

(2002). 

M. Baiesi and S. S. Manna, arXiv: cond-mat/0305054 

A. Capocci, G. Caldarelli, R. Marchetti^ and L. 

Pietronero, Phys. Rev. E. 64 R035105 (2001). 

D. A. Fell and A. Wagner, Nature Biotechnology 18, 

1121(2000). 

R. V. Sole, R. Pastor-Satorras, Handbook of Graphs and 
Networks ( Wiley- VCH, Berlin, 2003). 



B. Tadic, Physica 



